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ABSTRACT 



Context. Strong, globally-organized magnetic fields are found for a small fraction of O, B, and A stars. At the same time, many the- 
oretical and indirect observational studies suggested ubiquitous presence of weak localized magnetic fields at the surfaces of massive 
stars. However, no direct detections of such fields have been reported yet. 

Aims. We have carried out the first comprehensive theoretical investigation of the spectropolarimetric observational signatures of the 
structured magnetic fields. These calculations are applied to interpret null results of the recent magnetic surveys of massive stars. 
Methods. The intensity and circular polarization spectra of early-type stars are simulated using detailed polarized radiative transfer 
calculations with UTE model atmospheres. Similar to observational analyses, the mean Stokes / and V line profiles are obtained by 
applying a multi-line averaging technique. Diff'erent spectropolarimetric observables are examined for multiple realizations of ran- 
domly distributed radial magnetic field spots with diff'erent spatial scales. 

Results. We characterize the amplitude of the circular polarization profiles and the mean longitudinal magnetic field as a function of 
magnetic spot sizes. The dependence of these observables on the eff'ective temperature, projected rotational velocity, and inclination 
angle is also investigated. Using results of the recently completed Magnetism in Massive Stars (MiMeS) survey, we derive upper 
limits on the small-scale magnetic fields compatible with the MiMeS non-detections. 

Conclusions. According to our simulations, existing spectropolarimetric observations of sharp-lined massive stars rule out the pres- 
ence of the small-scale fields stronger than 50-250 G, depending on the typical spot sizes. For broad-lined stars, the observations 
constrain such fields to be below approximately 1 kG. 
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1. Introduction 

There is a growing observational and theoretical support for 
a magnetic dichotomy among intermediate-mass and massive 
stars. A small fraction of the late-B and A stars possesses 
strong organized magnetic fields at their surfaces (see reviews by 
iDmiati & Landstreet, 2009; Mathys, 2009; Kochukhov, 2011). 
These magnetic field topologies are approximately dipolar for 
. the majority of stars, often inclined relative to the stellar rota- 
' tional axis, and have strengths of up to ~ 30 kG. The global mag- 
netic fields of intermediate-mass stars are believed to be stable 
fossil remn ants of the magnetic flux acquire d at the stellar f orma- 
I tion phase (iBraithwaite & Nordlund[l2QQ6b . These fields exhibit 
. a lower threshold of ~ 300 G, below which the large-scale mag- 
■ netic topologies are probab ly not able to withstand the shearing 
by the diflTerential rotation (lAuriere et al.i I2007D . 

The incidence of fossil stellar magnetic fields is now 
understood to extend well into the mass range of early- 
B and s tars. T he Magnetism in Massi ve Stars (MiMeS ; 
'Wade etal.","2011a') survey and other studies (Wade et all l2006l: 
Hubrig et al., 2011) have identified several examples of O stars 
with globally-organized, kG-strength magnetic fields. These dis- 
coveries are particularly frequent among the so-c alled Of?p stars 
(iMartins et al.L 120101: IWade et al.L l2011bi l2012ah . which appear 
to represent a massive-star extension of the Ap/Bp phenomenon. 
Remarkably, magnetic fields existing on the surfa ce of some of 
these peculiar O stars exceed 20 kG jWade et al.L [2012cb , rival- 
ing the strongest fields found in much smaller Ap/Bp stars. 



In addition to the strong, large-scale magnetic fields found 
for a handful of massive stars there exists a substantial indirect 
observational evidence and credible theoretical predictions for 
the existence of localized, and probably weak, magnetic fields 
in the majority of massive stars. Much of the indirect support 
for such fields is coming from various poorly understood wind 
and surface variability phenomena. For example, nearly all O- 
type stars and many B st ars show unexplai ned cyclical vari- 
ability in th e ir win ds (Kap er et al .Lfl996: Fullerton et al. I [19961: 
iPrinja et all l2002h . The most prominent feature of this vari- 
ability is discrete absorption components (DACs), which are 
observed to accelerate through the UV-wind line profiles and 
maintain coherency over several stellar rotations. These features 
are probably caused by corotating interactin g regions (CIRs ) in 
the line-driven stellar wind, as proposed by iMuUanl (Il984h . In 
this model, a faster moving stream collides with a slow mov- 
ing one, leading to the formation of curved large scale structure 
in the stellar wind. The modeling by lCranmer & Qwockil(ll996h 
demonstrated that such wind structures give rise to DACs in the 
wind- sensitive UV line profiles. The velocity contrast required 
by this model arises from perturbations at the photospheric level, 
which may be connected to evolving small-scale magnetic spots. 

An example of the compelling observational evidence of 
the connection between DACs and a magnetic field was pro- 
vided by the long-term monito ring of the magnetic early-0 star 
0^ Ori C. IHenrichs et al.l (l2005b showed that 13 years of the C iv 
UV profile observations of this star, phase-folded with the 15.4 d 
rotational period, show a progressing absorption feature very 
similar to individual DACs observed in other early-type stars. 
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0^ Qri C is an oblique rotator with a sta ble dipolar magnetic field 
(iDonati et al.l.l2QQ2l: IWade et al.' '2006'). so the feature in the UV 
wind lines reoccurs exactly at the same rotational phase, asso- 
ciated with one of the magnetic poles. Using this well- studied 
0-star as an example, one can suspect magnetic field to be the 
main culprit of the surface perturbations leading to DACs in 
other massive stars. 

A considerable theoretical eflTort was devoted to investi- 
gating various alternatives to the fossil magnetism of early- 
type stars. Massive st ars generate stron g magnetic fields in 
their convective cores (iBrun et al.l l2005l) . but these fields do 
not reach the st ellar surface on the observable time scales 
(iMacGregor & Cassinelli . 2003). A dynamo process op erating 
in the diflTerentially-rotating radiative stellar interiors (ISpruitl 
l2QQ2h can also produce magne tic fields e merging on the 
stellar surface (Mullan & MacDonaldl l2005l) . The reality of 
this particular magnetic field generation mechanism, which 
has important impl i cation s for the evolution of r nassive stars 
(iMaeder & MevnetL l2005h and t heir descendants (iHeger et al.L 
120051) . was however disputed by'Zahn etal.' (2001). Finally, re- 
cent stellar structure calculations of hot massive stars suggested 
the presence of sub- surface convectiv ely unstable re gions, re- 
lated to the iron opacity peak ( Cantiello et al.L l2009l) . By anal- 
ogy with the turbulent magnet o-convection operati ng in other 
astrophysical environments (e.g.'Kapyla et al.V2008), these sub- 
surface convection zones (SCZ) can produce magnetic fields in 
the presence of shear and stellar rotation and provide an at- 
tractive framework for the explanation of other massive-star at- 
mospheric phenomena, includi ng excessi v e micr o- and macro- 
turbulent line broadening (Hunter et alj l2008l) and stochas- 
tic non-radial oscilla tions (Degroote et al. i I2010h . Preliminary 
MHD simulations by ICantiello et al.l J201 ll) confirmed the abil- 
ity of SCZ to generate localized magnetic fields visible at 
the stellar surface. The properties of such fields and associ- 
ated brightness spots were fur ther investigated analytically by 
ICantiello & Braithwait ^ (l201 ih . 

Simultaneously with these developments for massive stars, 
the discovery of sub-G magne tic fields in the brightest 
intermediate-mass stars Vega (Lign ieres et al.l l2009h and Sirius 
(iPetit et al.l. 1201 1) hinted at the hidden magnetism of nor- 
mal A and late-B stars. I n response to these findings, 
iBraithwaite & Cantiellol (1201 3l) developed a theoretical model of 
the dynamically-evolving "failed" fossil magnetic fields, which 
are expected to be much weaker and structured on smaller scales 
than the well-known stable fossil fields of the magnetic Ap/Bp 
stars. 

Despite an emerging interest in probing the small-scale mag- 
netism of early-type stars and exploring the consequences for the 
stellar physics, the current theoretical models have not reached 
the level of sophistication and detail necessary for making mean- 
ingful predictions about the strengths and topologies of the lo- 
calized magnetic fields in individual stars. Therefore, currently 
there exists a gap between observational constraints, e.g. detec- 
tions of the polarization signatures for Sirius and Vega and non- 
detections of the magnetic fields for most O- and B-type stars 
included in the MiMeS survey, and theoretical predictions. In 
particular, it is unknown what direct observational signatures 
of small-scale fields one can expect in the polarization spectra 
of the early-type stars and how the amplitude and variability of 
these features change with the characteristics of magnetic field 
and with the stellar parameters. 

Even the very ability of modern spectropolarimetric methods 
to detect localized magnetic features on the stellar surfaces is a 
matter of debate. On the one hand, a contribution of the small- 



scale fields to the disk-integrated polarization spectra and to the 
resulting magnetic observables might be strongly reduced by the 
cancellation of diflTerent field polarities. On the other hand, high- 
resolution studies of polarization signatures inside spectral line 
profiles partly overcomes this problem t hanks to the r otational 
Doppler resolution of the stellar surface (ISemeli Il989l) . The in- 
terplay of the cancellation vs. Doppler resolution and its depen- 
dence on the stellar rotation and sizes of magnetic features has 
not been systematically investigated. 

The goal of our study is to address several of these issues 
through a comprehensive and realistic numerical simulations of 
the circular polarization spectra corresponding to the small-scale 
magnetic fields of diflferent degree of complexity. In the follow- 
ing, we describe the methodology of the magnetic spectrum syn- 
thesis and multi-line analysis of the theoretical Stokes V spec- 
tra (Sect. |2]), present a summary of the simulation results for a 
wide range of stellar and magnetic field parameters (Sect. |3} and 
then use these data to constrain the strength of possible small- 
scale stellar magnetic fields by comparing our predictions with 
the outcome of recent magnetic field surveys of massive stars 
(Sect. a. 

2. Methods 

Two main spectropolarimetric observational techniques are cur- 
rently being used in the context of large-scale magnetic field sur- 
veys of early-type stars. The low-resolution circular spectropo- 
larimetry, using e.g. the FORSl/2 instrument at the ESO VLT, 
detects magnetic field by measuring the mean longitudinal mag- 
netic field in the wings of the hydrog en B aimer lines or in the un- 
resolved blends of metal l i nes (Bagnulo et al.Ll2QQ2l: IWade et all 
120071: iHubrig et al.L l201 ll: iBagnulo et al.1. l2012l)~ This technique 
is only sensitive to the line of sight magnetic field compo- 
nent integrated over the visible stellar hemisphere and is poorly 
suited to provide constraints on intermittent, small-scale mag- 
netic structures. Furthermore, flexures and calibration problems 
of the Cassegrain-mountedpolarimeters, such as FORSl/2, limit 
their usefulness for the search of very weak magnetic field (see 
discussion in Bagnulo et al., 2012). 

The second common magnetic field detection technique re- 
lies on obtaining high-resolution circular polarization spectra 
covering a wide wavelength range and applying a line- averaging 
technique to combine information from all suitable metal ab- 
sorption features. This approach is adopted by the MiMeS sur- 
vey for the analysis of circular polarization spectra obtained 
with ESPaDOnS, NARVAL and HARPSpol sp e ctropolarime- 
ters (iGrunhut et al.L 120091: lAlecian et all 120111: iNeiner et al.l 
l20T2h . The mean Stokes V line profiles obtained with the least- 
squares deconvolution (LSD) method (Donati et al.. 1997.) are 
particularly useful for the magnetic fi eld diagnostic since they 
behave similarly to real spectral lines ('K ochukhov et al.L 120101) 
and allow one to characterize and model spectrally-resolved 
Zeeman signatures. The LSD profiles are therefore sensitive to 
the field to pologies yielding v ery low or null mean longitudinal 
fields (^e.g. IShultz et al.l l2012h and have been successfully used 
to detect and map complex dynamo-generated magnetic fields of 
cool active stars (Donati et al., 2003; Kochukhov et al., 20131). 

The goal of our study is to characterize observational signa- 
tures of the magnetic fields which are at least as complex topo- 
logically as the surface magnetic structures found for cool ac- 
tive stars. Hence, the high-resolution spectropolarimetry com- 
bined with a multi-line analysis appears to be the most suitable 
methodology for our investigation. More specifically, we simu- 
late the LSD Stokes / and V profiles for magnetic field geome- 
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tries of varying degree of complexity and study these polariza- 
tion profiles as a function of the magnetic and stellar parame- 
ters. The following subsections describe computation of the lo- 
cal LSD line profiles, explain how we generate random magnetic 
field maps and perform surface integration and, finally, discuss 
derivation of different spectropolarimetric observables from the 
resulting disk-integrated polarization profiles. 

2.1. Local polarized line profiles 

The analysis of the LSD profiles of early-type stars presented 
in previous spectropolarimetric studies relied on a simplified 
analytical description of the intensity and circular polarization 
profile shapes. The most common assumption is to describe 
Stok es / as a Gaussian a nd Stokes V as its scaled derivative 
(e.g. iPetit & Wad3. l2Q12h . Although straightforward to imple- 
ment, this approach suffers from ambiguities in the choice of an 
analytical line profile model and its parameters, and is unable to 
properly describe the temperature dependence of the local line 
profile characteristics. Moreover, averaging over a small number 
of diagnostic lines available in the optical spectra of early-type 
stars does not necessarily yield a simple, symmetric Stokes / 
profile that could be described by an analytical function. Here 
we avoid these problems by numerically deriving the local LSD 
profiles from detailed polarized radiative transfer calculations of 
individual spectral lines. 

The local Stokes / and V line profile tables were prepared by 
first computing the full local spectra covering the 4000-8000 A 
wavelength range with the synmast polarized radiative transfer 
code and then applying the LSD line averaging procedure. A 
detailed description of the relevant computer codes and method- 
ology requir ed for this analys i s can be found in our previous 
publications (iKochukhov et al.Lr201Ql) . 

The spectrum synthesis calculations were performed for five 
ATLAs9 LTE model atmospheres (Kurucz, 1993) with Teff in the 
range from 10000 to 30000 K with a step of 5000 K. For all 
models we adopted the surface gravity of log ^ = 4.0 and micro- 
turbulent velocity of 2 kms"^. The atomic line dat a were ex- 
tracted from the VALD database (Kupka et al.Vl999'), using the 
1% line depth selection threshold and the same microturbulent 
velocity. The spectral regions containing the hydrogen and he- 
lium lines were subsequently excluded from the input line lists. 

The local Stokes IV parameters and the unpolarized contin- 
uum intensity /c were computed for a grid of 20 limb angles 6, 
distributed equidistantly with [i = cos 6, and for the two val- 
ues of local magnetic field strength, B = and 0.5 kG. For the 
radial magnetic field geometries considered below, the line of 
sight magnetic field component required for the polarized spec- 
trum synthesis is given by 5z = Bfi, independently of the dis- 
tribution of magnetic features over the stellar surface. For all 
spectrum synthesis calculations we employed a microtubulence 
of 2 kms"^ consistently with the choice of model atmospheres 
and the VALD line list extractions. 

On the next step, the local Stokes I(A,/u,B) and V(A,/u,B) 
spectra were normalized by the corresponding continuum inten- 
sity Ic(A,/j,) and convolved with the instrumental profile corre- 
sponding to the resolution of 7? = 65 000, appropriate for the 
ESPaDOnS and NARVAL spectropolarimeters. The resulting 
spectra were converted to the local LSD profiles, /lsdCv,//, 5) 
and VlsdCv,//, 5), by applying the least-squares deconvolution 
procedure to a subset of spectral lines deeper than 10% of 
the continuum. Simultaneously, the mean continuum intensity 
Ic(fi) was established by averaging Ic(A,fi) in the central wave- 
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Fig. 1. Normalized disk-center LSD Stokes / and V profiles for 
diff'erent eflTective temperatures and a line-of- sight magnetic field 
of 0.5 kG. The Stokes V spectra are oflfset vertically by 1.05 and 
expanded by a factor of 5. 

Table 1. Characteristics of the line lists used for synthetic LSD 
profile calculations. 
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lengths of the spectral lines included in the LSD masks. The 
LSD profiles were computed on a fine velocity grid with a step 
of 0.5 km s"^ covering the ±500 km s"^ velocity range to ensure 
an accurate surface integration for Ve sin / of up to 250 km s"^ . 

Table [T] summarizes statistics of the line lists used for the 
spectrum synthesis and LSD profile calculations. Depending on 
the temperature, we used 1500-2600 lines for the spectrum syn- 
thesis and up to ^ 400 lines for LSD. The mean wavelength Aq, 
eflTective Lande factor g^, and line depth do of the LSD masks 
are reported in the last three columns of Table [T] The same mean 
parameters were us ed to renormalize the LSD weights for each 
Teff as described by IKochukhov et af] (12010). 

Examples of the local disk-center LSD Stokes / and V pro- 
files for reff = 10000, 20000, and 30000 K are presented in 
Fig.[T] The strongest LSD profiles are obtained for low Tqq. The 
line width systematically increases towards higher temperatures. 
The local Stokes parameter line shapes are not exactly symmet- 
ric and contain weak features outside the main profiles due to 
blends unaccounted for in the LSD mask and due to limitations 
of the line averaging procedure itself. These LSD profile distor- 
tions are, of course, also present in observations but cannot be 
described by an analytical line profile model. 

2.2. Disk-integrated spectra for random magnetic fields 

Many diflferent mathematical models can be used to construct 
small-scale magnetic field distributions on the stellar surfaces. 
Here we employ a simple formulation, requiring few arbitrary 
parameters but still maintaining sufliicient physical realism. In 
the following we assume that the stellar magnetic field topol- 
ogy is given by a superposition of randomly distributed circular 
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Fig. 2. Examples of random distributions of radial magnetic field spots with filling factor /sp = 0.5 and rgp = 10° {left panel) and 2° 
{right panel). In each case the star is shown at nine rotational phases and inclination angle / = 60°. The white and black regions 
correspond to the spots with positive and negative field polarity, respectively. The grey areas are field-free. 



spots. The field inside the spots is purely radial, by analogy with 
sunspots, and is characterized by a single value of the magnetic 
field modulus. Outside the spots the photosphere is assumed to 
be non-magnetic. An equal number of spots with the outward- 
and inward-directed radial field is used, thus maintaining the 
zero net magnetic flux through the stellar surface. The average 
properties of such magnetic field distributions are defined by 
three parameters: the spot angular radius rsp, the field strength 
^sp inside the spots, and the fraction of the stellar surface /sp 
covered by the magnetic field (a filling factor). 

Once Tsp is set, the area of each spot can be computed as 



^sp 



27lR^{\ - COSTsp) 



(1) 



and the total number of spots for a given /p can be estimated by 
rounding 



^sp = 47TR^fsp/a,^ 



(2) 



to the nearest even integer. Then, the latitudes x ^ [-7r/2,7r/2] 
and longitudes (f e [0, 27i] of the spot centers are obtained from 



X = 7r/2 - arccos {2x - 1) 
(f = 2nx\ 



(3) 



where x and x' represent two different sequences of uniformly 
distributed random numbers in the range between and 1 . This 
algorithm of the random spot position generation may produce 
overlapping spots. In practice, we allowed the spots to overlap 
by no more than 10% of their area, rejecting spot coordinates 
too close to the already existing spots. 

The local Stokes / line profiles in Fig. [T] have a full- width at 
half maximum of = 8-14 km s"^ The corresponding approxi- 
mate theoretical spatial resolution limit can be estimated as 



71 Sv 



^ min — . 



4 Ve sin / 



(4) 



Using = 10 kms \ we can simplify this expression to rmin = 
900/ Ve sin/. This indicates that an isolated surface feature with 



Tsp = 4.5° is fully resolved for Ve sin / = 200 kms"^. Using this 
estimate as an guideline, we considered the spot sizes 2, 5, 10, 
20, and 40°. An intermediate magnetic filling factor /p = 0.5 
was adopted for the majority of calculations. The corresponding 
total spot numbers varied from 4 (for rsp = 40°) to 1642 (for 
rsp = 2°). 

Numerically, magnetic maps were evaluated on the spheri- 
cal surface gr id containing 15130 zones of ro ughly equal areas 
(see Fig. 5 in lPiskunov & Kochukhovl [20021) . This fine surface 
grid was required to properly resolve the spots with the smallest 
spatial scales. An example of the random magnetic spot distri- 
butions computed for rsp = 10 and 2° is illustrated in Fig.[2l 

After establishing a surface magnetic field map, we com- 
puted the disk-integrated Stokes / and V LSD profiles, (/lsd) and 
(^lsd), by appropriate interpolation and weighted summation of 
the Doppler- shifted local LSD Stokes profiles for a specific in- 
clination angle / and projected rotational velocity Ve sin /. The 
line profiles were simulated for 20 equidistant rotational phases, 
using a velocity grid with 1.8 kms"^ step size. All calculations 
presented below were performed for the Ve sin / values of 10, 20, 
50, 100, and 200 km s'^. 

An example of the phase variation of the circular polarization 
profiles for the random magnetic field distribution with rsp = 10° 
and /sp = 0.5 (the same as shown in the left panel of Fig. O is 
given in Fig. [3] and IH These figures illustrate the line profiles 
computed for identical surface maps and the same set of view- 
ing geometries, but diflferent projected rotational velocities. As 
expected, a higher projected rotational velocity (Fig. |4]l reveals 
much more complex circular polarization profile shapes thanks 
to the Doppler resolution of the stellar surface. At the same time, 
the amplitude of the Stokes V signatures is also considerably re- 
duced compared to low the Ve sin / case (Fig. [3]). 

All our calculations were performed for the 0.5 kG field 
strength inside magnetic spots. Since here we deal with the mag- 
netic line formation in the weak-field regime, it can be assumed 
that the Stokes V profiles and all resulting spectropolarimet- 
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Fig. 3. Disk-integrated normalized LSD Stokes / and V profiles 
for Ve sin / = 10 km s"^ and a random distribution magnetic spots 
with Tsp = 10°. Other parameters are T^^ = 20000 K, / = 60°, and 
magnetic filling factor /sp = 0.5. The spectra corresponding to 
diflTerent rotational phases are oflTset vertically. 



ric observables scale linearly with the magnetic field strength. 
Consequently, after verifying this behavior for the local Stokes 
V LSD profiles, we did not vary the field strength inside mag- 
netic spots in any of the numerical experiments presented in this 
paper. 

2.3. Magnetic observables 

The main goal of our study is to assess the capability of re- 
cent high-resolution spectropolarimetric surveys of early-type 
stars to detect intermittent magnetic fields on their surfaces. 
Consequently, we consider the amplitude of the Stokes V sig- 
natures resulting from our numerical line profiles simulations 
as the primary magnetic field observable. For each set of cir- 
cular polarization profiles we calculated the mean peak-to-peak 
Stokes V amplitude Vmean by averaging the instantaneous am- 
plitude over 20 rotational phases. We also determined the max- 
imum amplitude Vmax by finding the rotational phase with the 
largest circular polarization signature. To obtain statistically 
meaningful results with our random magnetic spot models, we 
averaged both Vmean and Vmax over 100 diflTerent random realiza- 
tions of the magnetic field distributions for the same rgp and /sp. 
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Fig. 4. LSD Stokes / and V profiles computed with the 
same magnetic map and setup as in Fig. [S] except for 
Ve sin/= 100 kms"^. 



These calculations were repeated for a grid of five r^p and five 
Ve sin / values (see above), yielding 50000 individual synthetic 
line profiles for each Tq^, /p, / combination. 

In addition to considering the amplitude of the Stokes V sig- 
natures, we introduced another parameter characterizing com- 
plexity of the circular polarization profiles. This (integer) num- 
ber was established by counting how many times a given Stokes 
V profile intersects the zero line. This parameter also roughly 
corresponds to the number of lobes in the Stokes V profile mi- 
nus one. The final mean reported below for each rsp-Ve sin / 
combination was obtained by computing the median over 20 ro- 
tational phases and 100 random map realizations. We found that 
this simple diagnostic exhibits a smooth variation with stellar 
parameters and readily distinguishes complex and simple polar- 
ization signatures. For example, for the magnetic spot distribu- 
tion with Tsp = 10° illustrated in Fig. |2]^z gradually increases 
from 2 for Ve sin/= 10 kms"^ (line profiles in Fig. [S]) to 9 for 
Ve sin / = 100 km s"^ (line profiles in Fig. |4]). 

For completeness, we also determined the mean longitudinal 
magnetic field (5z) from each pair of the disk-integrated Stokes 
/ and V profiles. This was accomplished using the equation 



. r V<yLSD)dv 

<5z) = -7.145 X 10^^^^-^ , 

^0^0 J(l -</LSD))dv 



(5) 
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Fig. 5. Mean longitudinal magnetic field and the LSD Stokes V profile amplitude as a function of the magnetic spot radius and pro- 
jected rotational velocity. These calculations were carried out for Tqq = 20 000 K and magnetic filling factor /sp = 0.5. The top panels 
show maximum K^z)! for a sample of 20 rotational phases (top left) and K^z)! averaged over this set of rotational phases (top right). 
The bottom panels illustrate the corresponding maximum (bottom left) and mean (bottom right) circular polarization amplitude. 
Each point in this plot was computed by averaging over 100 diflferent random realizations of the magnetic spot distributions. 



where Aq in A and the resulting (5z) is in G. The mean wave- 
length Aq and eflfective Lande factor g^ are the same param- 
eters as applied above for normalizing the LSD line weights 
(see Table [TJ. Similar to the analysis of the Stokes V ampli- 
tude, the maximum and mean unsigned longitudinal magnetic 
fields, K5z)lmax and K5z)lmean, wcrc first obtained from a set of 
theoretical spectra computed for 20 rotational phases and then 
averaged over 100 diflferent random realizations of the magnetic 
field maps with the same r^p and fp. 

3. Results 

Our main set of the line profile calculations for diflTerent rgp and 
VeSin/ was generated assuming Teff = 20000 K, / = 60°, and 



/sp = 0.5. Here we present results obtained with this combina- 
tion of parameters and then assess the impact of changing the 
eflTective temperature, inclination angle, and the filling factor. In 
all cases calculations were performed for B^p = 500 G. 

The mean longitudinal magnetic field and the amplitude of 
the Stokes V profiles are plotted as a function of the magnetic 
spot size and projected rotational velocity in Fig. [5] We illus- 
trate the maximum K^z)! and the largest polarization amplitude 
that might be expected for the observations with a dense phase 
coverage (left column) and the phase-averaged values of these 
observables, representative of snapshot observations (right col- 
umn). 

As one can see from the upper panels of Fig. [5] the maximum 
longitudinal magnetic field changes from ~ 10 G for r^p = 2° to 
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Fig. 6. Number of intersections of the zero line in the Stokes V 
profiles as a function of the magnetic spot radius and projected 
rotational velocity. Each point is a median computed over 20 
rotational phases and 100 dififerent random realizations of mag- 
netic spot distributions. 



150-180 G for rgp = 40°. The corresponding K5z)lmean is approx- 
imately two times smaller. There is a small variation of K^z)! 
with Ve sin /, which is not expected for the purely geometrical 
definition of this magnetic parameter. This Ve sin / dependence is 
real and is explained by the presence of weak structures in the 
local LSD Stokes spectra outside main line profiles (see Fig.[T]). 
These structures contribute diflTerently to the disk-integrated pro- 
files depending on the range of rotational Doppler shifts, result- 
ing in slightly diflTerent first-order moments of (Vlsd) for diflfer- 
ent Ve sin /. We have verified that the Ve sin / dependence of K^z)! 
disappears if we substitute numerical local LSD profile tables 
with (less realistic) smooth analytical profiles. 

The lower panels in Fig. [5] illustrate the mean and maximum 
Stokes V amplitude as a function of r^p and Ve sin /. This figure 
reveals a non-trivial interplay between the Doppler resolution 
and cancellation of the opposite spot polarities. For sharp -lined 
stars there is a strong dependence of the polarization amplitude 
on Tsp, with the largest spots yielding about 10 times higher am- 
plitude of the Zeeman polarization signatures. As soon as Ve sin / 
exceeds -100 km s"\ the amplitude of polarization profiles be- 
comes almost independent of r^p. The typical peak-to-peak line 
polarization amplitude obtained in our numerical experiments 
lies in the range from ~ 10"^ to ~ 10"^. 

According to our results, the low-Ve sin / stars show a higher 
amplitude polarization features independently of the typical 
magnetic spot size. In other words, an increase of the polariza- 
tion amplitude due to the rotational Doppler resolution does not 
compensate the reduction of the signal due to weakening of spec- 
tral lines. This suggests that massive sharp-lined stars (intrinsi- 
cally slow rotators or fast rotators visible pole-on) are preferred 
targets for searching for signatures of complex magnetic fields. 
However, it should be remembered that an observational assess- 
ment of the typical sizes of magnetic features is only possible 
for fast rotators. To this end. Fig. [6] presents an analysis of the 
Stokes V profile complexity using the number of zero line in- 
tersections, Hz, for the same range of Ve sin / and r^p as shown 
in Fig. \5\ Clearly, a high Ve sin / is beneficial for ascertaining 
the field complexity and resolving the surface structures at the 
smallest spatial scales. At the same time, a moderate Ve sin / of 
~ 50 km s"^ is already suflftcient for distinguishing, with a single 



high-quality spectropolarimetric observation, small-scale mag- 
netic fields from the global fossil magnetic topologies. 

All calculations presented here assumed 500 G field strength 
inside magnetic spots. Our results can be easily scaled to an ar- 
bitrary starspot field strength by multiplying the magnetic ob- 
servables by B^p/SOO (for the field strength measured in G). 
For the logarithmic plots shown in Fig. [5] this re-scaling corre- 
sponds to a vertical translation of all curves by lg(5sp/500). For 
example, decreasing the spot field strength to 100 G will result 
in the downward shift of the Stokes V profile amplitude curves 
by a factor of 0.7. The morphology of the circular polarization 
profiles and Fig. [6] will remain unchanged. 

In the following we describe the impact of variation of Tqq, 
/, and /sp with respect to the default values of these parameters 
(Jeff = 20000 K, / = 60°, and f,p = 0.5). First, repeating the full 
set of calculations for / = 30° and 90°, we found no dependence 
of any of the considered Stokes V profile characteristics on the 
stellar inclination. Statistically, the observational signatures of 
localized magnetic fields turn out to be completely insensitive to 
this parameter. 

The impact of using the local line profile tables correspond- 
ing to diflTerent Tq^ can be qualitatively assessed with the help of 
the local Stokes / and V profiles shown in Fig. [T] Compared to 
the full Teff range of 10000-30000 K, T^^ = 20000 K adopted 
for our main set of calculations represents the least favorable 
case due to the smallest Stokes / line depth. Unsurprisingly, us- 
ing Teff = 10000 K increases the (Vlsd) amplitude by 30% for 
rapidly rotating stars and up to 70% for sharp-lined stars. The 
change associated with increasing Tq^ from 20 000 to 30 000 K 
is less significant. In this case, an increase of the line depth is 
counteracted by an increase of the local profile width. As a re- 
sult, we obtained ^ 20% Stokes V signal enhancement for sharp- 
lined stars and no change or even a few per cent decrease of the 
polarization signal for rapid rotators. 

Finally, we investigated the impact of reducing the filling 
factor of the magnetic spots from 0.5 to 0.2. As expected, this 
large decrease of the amount of unsigned magnetic flux on the 
stellar surface yields some reduction of both the Stokes V am- 
plitude and the mean longitudinal magnetic field inferred from 
the circular polarization profiles. However, this reduction is not 
linear with /sp because a smaller number of magnetic spots on 
the visible hemisphere of the star also leads to a less eflftcient 
cancellation of the opposite magnetic field polarities. On aver- 
age, we found about 50% decrease of Vmax, Vmean and as well as 
K^z)lmax and K5z)lmean for a factor of 2.5 reduction of the surface 
area covered with magnetic spots. 

4. Discussion 

In the previous section we calculated the amplitudes of the mean 
longitudinal magnetic field and the circular polarization inside 
spectral lines for a fixed magnetic field strength and for diflTerent 
magnetic spot sizes. These predictions cannot be directly com- 
pared to any observations because no detection of small-scale 
magnetic fields has been reported so far for early-B and O stars. 
Nevertheless, it is of interest to interpret our results in terms of 
the maximum strength of such fields compatible with the most 
sensitive of the existing observational non-detections. 

The most comprehensive search of the magnetic fields in 
massive stars was carried out by the MiMeS survey using 
the ESPaDOnS, NARVAL, and HARPSpol spectropolarimeters. 
Several discoveries of the strong, large-scale magnetic fields 
have been made from these observations (e. g. iGrunhut et al.l 
120091: lAlecian et al. 1 . l201ll: IWade et al.L l2012cl) . In addition, the 
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Fig. 7. Upper limits for the small-scale magnetic field intensity 
corresponding to the null results of the MiMeS survey. Upper 
panel: B^^ corresponding to K^z)! = 30 G (solid lines) and 75 G 
(dashed lines). Lower panel: ^sp required for the magnetic field 
detection with the LSD Stokes V profiles at a confidence level 
of > 99.9% for at least half of the rotational phases. The dotted 
horizontal lines show ^sp = 50, 100, 200, 500, and 1000 G. 



MiMeS survey component (SC) include over 300 O and B stars 
for which no eviden ce of magnetic field was found. According to 
IWade et al.l (l2012bh . the resulting incidence of global magnetic 
fields is about 6% for both O and B stars. 

A summary of the magnetic field non-detections for the 
MiMeS SC stars was kindly communicated to us by J. Grunhut. 
The median uncertainty of the {B^} measurements is 25 G for all 
SC stars. Somewhat lower error bars of ~ 10 G were achieved 
for about two dozen stars with VeSin/< 100 kms"^ The me- 
dian signal-to-noise ratio of the LSD Stokes V profiles is 23 000. 
Using these observational constraints, we inferred the upper lim- 
its for the magnetic field strength using our line profile and K^z)! 
predictions for T^^ = 20000 K, ^ = 0.5, and / = 60°. 



The upper panel of Fig. [7] shows the minimum ^sp necessary 
for a 3cr detection of the mean longitudinal magnetic field for 
(T = 25 G and 10 G. Evidently, (5z) is not a particularly useful 
diagnostic when it comes to small-scale magnetic field topolo- 
gies. The field detection using (5z) requires ^sp ^ 500 G for a 
subgroup of the MiMeS SC targets with the most accurate longi- 
tudinal magnetic field measurements. This limit is B^p ^ 1500 G 
for the entire SC sample. 

Next, we assessed the ability of the MiMeS survey to detect 
small-scale magnetic fields with the LSD line profile analysis. 
We conducted a series of Monte-Carlo simulations by adding 
a random noise with diflTerent amplitude to the theoretical LSD 
Stokes V profiles corresponding to r^p = 10° and estimating the 
probabilit y of polarizat ion signal detection using the chi-square 
statistics (iDonati et~aD. 1199 2). These calculations showed that, 
independently of the Stokes V profile morphology, the circular 
polarization signatures can be detected at a confidence level of 
> 99.9% for at least half of the rotational phases if the noise 
is at least five times smaller than the peak-to-peak amplitude 
of the Stokes V profile. Combining this estimate with a typi- 
cal precision of the Stokes V LSD profiles obtained by MiMeS, 
we could infer the maximum possible strength of the magnetic 
field inside small-scale spots. This ^sp limit is illustrated in the 
lower panel of Fig. [7] as a function of r^^ and Ve sin /. This fig- 
ure confirms that the LSD Stokes V profile analysis is consider- 
ably more sensitive to the small-scale fields compared to (^z)- 
We find that for sharp-lined stars the detection of ^sp ~ 50 G 
field is possible for r^p > 10°. The limiting field strength in- 
creases to 200-300 G for the smallest spatial scales we have 
considered (r^p = 2°). Fast-rotating stars show a less pronounced 
dependence of the maximum ^sp on r^p. Typically, one should 
be able to detect the signatures of 300-500 G magnetic fields 
for any r^p at Ve sin/ ==100 kms"^ and 800-1200 G fields at 
Ve sin / = 200 km s"^ Interestingly, for fast rotating stars the limit 
of 800 G corresponds to the field with rsp < 5° and it is compar- 
atively more difl&cult to detect larger- scale magnetic field struc- 
tures. 

How do these constraints compare with theoretical predic- 
tions? The only quantitative estimate of the strength of localized 
magnetic fiel ds emerging on the surfaces o f massive stars was 
presented by Cantiello & Braithwaitel (l201ll hereafter CB 1 1 ) in 
the context of their study of the effects of turbulence in the iron 
SCZ. Supposing that the dynamo action in the convection zone 
is in equipartition with the kinetic energy of the gas, the authors 
found magnetic fields of up to ~ 2 kG at the base of the sur- 
face radiative zone. Further assuming that the magnetic fields are 
brought up to the surface by the magnetic buoyancy, CB 1 1 ob- 
tained the magnetic field strength of the order of 20-150 G in the 
photospheric layers, depending on the magnetic flux tube geom- 
etry and stellar parameters. CBl 1 also did not rule out the pho- 
tospheric fields reaching the equipartition limit of about 300 G. 

Little is known about the typical sizes of the surface mag- 
netic field structures generated by the SCZ. CB 1 1 argued that the 
lower limit of the magnetic starspot size is given by the pressure 
scale height, which is about 0.2 Rq for a massive star with R^, = 
IORq. This corresponds to rsp ^ 1°, i.e. close to the lower spa- 
tial limit investigated in our paper. At the same time, magneto- 
hydrodynamic simulations demonstrate formation of magnetic 
structures o n the scales significantly larger than the scale of 
convection (iKapvla et al.1. l2008l : ICantiello et"aD . l201ll) . For the 
Sun, the observed active regions also have sizes significantly 
exceeding the pressure scale heig ht (e.g. Zhang eLal., 2010) . 
Inhomogeneous wind models by dCranmer & Owockii 1 19961) 
considered surface spots with rsp ^ 10°, hence the magnetic 
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concentrations with coherent field orientation should be simi- 
larly large for the magnetic interpretation of CIRs and DACs. 
Taking into account these arguments, we suggest that r^p = 5- 
10° might be a more realistic estimate of the typical size of the 
hypothetical surface magnetic field regions. 

A remarkable result of our analysis is that the high-resolution 
Stokes V observations of massive sharp-lined stars have already 
reached the precision necessary for detecting r^^ = 5-10° spots, 
filled with ^sp = 50-100 G radial magnetic field and occupy- 
ing half of the stellar surface. However, no discovery of a corre- 
sponding population of narrow-line magnetic massive stars have 
been reported by MiMeS. This may indicate that the theoretical 
predictions by CBU overestimate the surface field strength or 
that the fraction of the stellar surface covered by magnetic re- 
gions is much smaller than f^^ = 0.5 adopted in our simulations. 
A more firm conclusion emerging from our work is the absence 
of equipartition fields with B^p ~ 300 G on the surfaces of mas- 
sive stars. 

The tentative comparison of the overall MiMeS survey re- 
sults with the generic theoretical predictions for a large range of 
stellar parameters can be put on more solid footing by consider- 
ing fundamental parameters, observational constraints, and the- 
oretical predictions for individual stars. On the one hand, such 
analysis would take into account the individual line lists, Ve sin /, 
and circular polarization detection limits and, on the other hand, 
would allow choosing an appropriate maximum theoretical mag- 
netic field strength according to the stellar luminosity and tem- 
perature (see Fig. 3 in CBll). In the future, it will be straight- 
forward to perform this analysis of individual massive stars with 
the tools and methodology developed in our paper. 

Acknowledgements. We thank Dr. Jason Grunhut for providing a summary of 
the MiMeS null results ahead of publication. OK is a Royal Swedish Academy 
of Sciences Research Fellow, supported by the grants from Knut and Alice 
Wallenberg Foundation and Swedish Research Council. NS is supported by the 
Saint-Petersburg State University through a research grant 6.38.73.2011. 



Kochukhov, O. 2011, in lAU Symposium, Vol. 273, lAU Symposium, ed. 

D. Prasad Choudhary & K. G. Strassmeier, 249-255 
Kochukhov, O., Makaganiuk, V., & Piskunov, N. 2010, A&A, 524, A5 
Kochukhov, O., Mantere, M. J., Hackman, T., & Hyin, I. 2013, A&A, 550, A84 
Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C, & Weiss, W. W. 

1999, A&AS, 138, 119 
Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s 

grid. Kurucz CD-ROM No. 13. Cambridge, Mass.: Smithsonian 

Astrophysical Observatory. 
Lignieres, F., Petit, R, Bohm, T., & Auriere, M. 2009, A&A, 500, L41 
MacGregor, K. B. & CassineUi, J. R 2003, ApJ, 586, 480 
Maeder, A. & Meynet, G. 2005, A&A, 440, 1041 

Martins, F., Donati, J., Marcolino, W. L. F., et al. 2010, MNRAS, 407, 1423 
Mathys, G. 2009, in Astronomical Society of the Pacific Conference Series, Vol. 

405, Solar Polarization 5: In Honor of Jan Stenflo, ed. S. V. Berdyugina, K. N. 

Nagendra, & R. Ramelh, 473 
MuUan, D. J. 1984, ApJ, 283, 303 

MuUan, D. J. & MacDonald, J. 2005, MNRAS, 356, 1139 

Neiner, C, Grunhut, J. H., Petit, V, et al. 2012, MNRAS, 426, 2738 

Petit, R, Lignieres, F., Auriere, M., et al. 2011, A&A, 532, L13 

Petit, V. & Wade, G. A. 2012, MNRAS, 420, 773 

Piskunov, N. & Kochukhov, O. 2002, A&A, 381, 736 

Prinja, R. K., Massa, D., & FuUerton, A. W. 2002, A&A, 388, 587 

Semel, M. 1989, A&A, 225, 456 

Shultz, M., Wade, G. A., Grunhut, J., et al. 2012, ApJ, 750, 2 
Spruit, H. C. 2002, A&A, 381, 923 

Wade, G. A., Alecian, E., Bohlender, D. A., et al. 2011a, in lAU Symposium, 
Vol. 272, lAU Symposium, ed. C. Neiner, G. Wade, G. Meynet, & G. Peters, 
118-123 

Wade, G. A., Bagnulo, S., Drouin, D., Landstreet, J. D., & Monin, D. 2007, 

MNRAS, 376, 1145 
Wade, G. A., FuUerton, A. W., Donati, J.-F., et al. 2006, A&A, 451, 195 
Wade, G. A., Grunhut, J., Grafener, G., et al. 2012a, MNRAS, 419, 2459 
Wade, G. A., Grunhut, J. H., & MiMeS Collaboration. 2012b, in Astronomical 

Society of the Pacific Conference Series, Vol. 464, Circumstellar Dynamics 

at High Resolution, ed. A. C. Carciofi & T. Rivinius, 405 
Wade, G. A., Howarth, I. D., Townsend, R. H. D., et al. 2011b, MNRAS, 416, 

3160 

Wade, G. A., Maiz Apellaniz, J., Martins, F., et al. 2012c, MNRAS, 425, 1278 
Zahn, J.-R, Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 
Zhang, J., Wang, Y., & Liu, Y. 2010, ApJ, 723, 1006 



References 

Alecian, E., Kochukhov, O., Neiner, C, et al. 2011, A&A, 536, L6 
Auriere, M., Wade, G. A., Silvester, J., et al. 2007, A&A, 475, 1053 
Bagnulo, S., Landstreet, J. D., Fossati, L., & Kochukhov, O. 2012, A&A, 538, 
A129 

Bagnulo, S., Szeifert, T., Wade, G. A., Landstreet, J. D., & Mathys, G. 2002, 

A&A, 389, 191 
Braithwaite, J. & Cantiello, M. 2013, MNRAS, 428, 2789 
Braithwaite, J. & Nordlund, A. 2006, A&A, 450, 1077 
Brun, A. S., Browning, M. K., & Toomre, J. 2005, ApJ, 629, 461 
Cantiello, M. & Braithwaite, J. 2011, A&A, 534, A140 

Cantiello, M., Braithwaite, J., Brandenburg, A., et al. 2011, in lAU Symposium, 
Vol. 272, lAU Symposium, ed. C. Neiner, G. Wade, G. Meynet, & G. Peters, 
32-37 

Cantiello, M., Langer, N., Brott, I., et al. 2009, A&A, 499, 279 

Cranmer, S. R. & Owocki, S. R 1996, ApJ, 462, 469 

Degroote, R, Briquet, M., Auvergne, M., et al. 2010, A&A, 519, A38 

Donati, J.-F., Babel, J., Harries, T. J., et al. 2002, MNRAS, 333, 55 

Donati, J.-F., Cameron, A. C, Semel, M., et al. 2003, MNRAS, 345, 1145 

Donati, J.-F. & Landstreet, J. D. 2009, ARA&A, 47, 333 

Donati, J.-F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, 

MNRAS, 291,658 
Donati, J.-R, Semel, M., & Rees, D. E. 1992, A&A, 265, 669 
FuUerton, A. W., Gies, D. R., & Bolton, C. T. 1996, ApJS, 103, 475 
Grunhut, J. H., Wade, G. A., Marcolino, W. L. R, et al. 2009, MNRAS, 400, L94 
Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 
Henrichs, H. R, Schnerr, R. S., & ten Kulve, E. 2005, in Astronomical Society of 

the Pacific Conference Series, Vol. 337, The Nature and Evolution of Disks 

Around Hot Stars, ed. R. Ignace & K. G. Gayley, 114 
Hubrig, S., SchoUer, M., Kharchenko, N. V, et al. 2011, A&A, 528, A151 
Hunter, I., Lennon, D. J., Dufton, R L., et al. 2008, A&A, 479, 541 
Kaper, L., Henrichs, H. R, Nichols, J. S., et al. 1996, A&AS, 116, 257 
Kapyla, R J., Korpi, M. J., & Brandenburg, A. 2008, A&A, 491, 353 



9 



